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The exploration of topological states is of significant fundamental and practical importance in contemporary 
condensed matter physics, for which the extension to two-dimensional (2D) organometallic systems is particu¬ 
larly attractive. Using first-principles calculations, we show that a 2D hexagonal triphenyl-lead lattice composed 
of only main group elements is susceptible to a magnetic instability, characterized by a considerably more stable 
antiferromagnetic (AFM) insulating state rather than the topologically nontrivial quantum spin Hall state pro¬ 
posed recently. Even though this AFM phase is topologically trivial, it possesses an intricate emergent degree of 
freedom, defined by the product of spin and valley indices, leading to Berry curvature-induced spin and valley 
currents under electron or hole doping. Furthermore, such a trivial band insulator can be tuned into a topo¬ 
logically nontrivial matter by the application of an out-of-plane electric field, which destroys the AFM order, 
favoring instead ferrimagnetic spin ordering and a quantum anomalous Hall state with a non-zero topological 
invariant. These findings further enrich our understanding of 2D hexagonal organometallic lattices for potential 
applications in spintronics and valleytronics. 

PACS numbers: 73.43.-f, 72.20.-i, 72.80.Le, 81.05.Fb 


Since the discovery of the quantum Hall effect |[ll], the con¬ 
cept of topological order has become a subject of major in¬ 
terest in contemporary condensed matter physics 12]. Kane 
and Mele revealed 12] that graphene, a two-dimensional (2D) 
honeycomb lattice of carbon atoms 10], opens a gap by spin- 
orbit coupling (SOC) and realizes a topologically nontrivial 
state called the quantum spin Hall (QSH) state, giving rise to 
a quantized response of a transverse spin current to an applied 
in-plane electric field. This discovery triggers a huge amount 
of activities in exploring 2D or 3D topological insulators (TIs) 
that possess robust helical conducting edge or surface states 
on the boundary of bulk insulators IdI-HQI] . This peculiar heli¬ 
cal state is topologically protected from elastic backscattering 
by time-reversal symmetry and hence offers fascinat¬ 

ing playgrounds for applications in spintronics and quantum 
computation devices Bill] . 


Analogous to the quantum control of the spin (HQ, other 
binary quantum degrees of freedom such as two inequivalent 
valleys (K and K' points) and two inequivalent lattice sites 
(A and B sublattices) in 2D hexagonal lattices, as shown in 
Fig. 1, have been exploited to bring the emergence of the 
so-called valleytronics iBlIl and pseudospintronics iIB- 
[H. As the inversion symmetry is broken in 2D hexagonal 
structures, Bloch wavefunctions exhibit opposite Berry cur¬ 
vatures at K and K' to host a quantum valley Hall (QVH) 
effect (il S m characterized by a valley Chern number 
Cv = Ck — Ck'- Recently, various schemes have been pro¬ 
posed to generate valley currents in graphene iQIIllIl and 
2D transition metal dichalcogenides 122, by using their 


unique edge modes |[l4l] . defect lines ll^ . and strain ll^ . 
Since a strong SOC gives rise to the interplay between spin 
and valley and lattice pseudospins O, it is interesting and 
challenging to examine how it can rearrange the band struc¬ 
ture and further change the topological property using exter¬ 
nal perturbations, e.g., electric field. Indeed, there have been 
several recent proposals on the topological phase transitions 
in 2D hexagonal structures, induced by an out-of-plane exter¬ 
nal electric field ll22J^ . However, the realization of these 
exciting possibilities in 2D materials is still in its infancy. 



FIG. 1: (Color online) (a) Top and side views of the structure of 
the 2D hexagonal triphenyl-lead lattice. The dotted line indicates the 
unit cell. The large, medium, and small circles represent Pb, C, and H 
atoms, respectively. Two inequivalent Pb lattice sites are denoted as 
A and B sublattices. The Brillouin zone is displayed in (b), together 
with the high symmetric ^-points F, M, K, and K'. 

Here we focus on a prototypical example of 2D hexagonal 
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organometallic systems, triphenyl-lead (TL) lattice ll^ . As 
shown in Fig. 1, Pb atoms form a honeycomb lattice with a 
threefold rotational symmetry by sharing phenyl bridges with 
three neighbors, and the A and B sublattices of Pb atoms ex¬ 
hibit slight puckering. Considering that the Pb atom has the 
6s^6p^ valence-electron configuration, one expects that the 
TL lattice would realize a QSH state, similar to silicene that 
is a puckered honeycomb lattice of Si S. Indeed, in their 
pioneering density-functional theory (DFT) studies of topo¬ 
logical phases in organic systems 1291 Isij], Wang, Liu, and 
Liu demonstrated that the TL lattice with a nonmagnetic con¬ 
figuration would be a QSH insulator. 

In this Letter, we perform comprehensive DFT calculations 
for the TL lattice to investigate its true ground state by includ¬ 
ing a quantum anomalous Hall (QAH) state as well as a band 
insulator with an antiferromagnetic (AFM) spin ordering. We 
find that the latter AFM insulating state is more stable than 
the QSH and QAH states. Despite its topologically trivial 
feature, this AFM ground state possesses a novel electronic 
degree of freedom defined by the product of spin and valley 
indices, leading to Berry curvature-induced spin and valley 
currents under carrier doping or polarized optical pumping. 
We also reveal that the application of vertical gating induces 
a topological phase transition from the AFM insulating state 
by closing one of the valley gap. As the gap reopens, the sys¬ 
tem enters into a ferrimagnetic QAH insulator with a Chern 
number of one. Because the AFM and QAH states with fas¬ 
cinating electronic properties are connected by a single ex¬ 
perimental knob and furthermore 2D organometallic systems 
are chemically flexible in terms of ligands and metal centers, 
the generic consequences of the present findings have much to 
offer for the design and development of spintronics and val¬ 
ley tronics applications. 

We begin to study the nonmagnetic (NM) state of the TL 
lattice using the PBE calculation without the inclusion of 
soc iH. The optimized structure is displayed in Fig. 1. 
The calculated distance Jpb-Pb between two neighboring Pb 
atoms and their height difference Ah are 7.452 and 2.932 A, 
respectively. Here, each phenyl ring rotates along the Pb-Pb 
axis by ^27° relative to the horizontal plane of an isolated 
TL molecule. The PBE band structure shows Dirac cones at 
the K and K' points, with a Dirac point at the Fermi level Ep 
(see Fig. lS(a) of the Supplemental Material Issl]). On the 
other hand, when SOC is included in the PBE calculation, the 
band gap Eg at K and K' is opened as high as ^ 11 meV Isee 
Fig. 2(a)], in good agreement with that (8.6 meV) of a pre¬ 
vious DFT calculation (2^. By computing the parity eigen¬ 
values of Bloch wavefunctions at the time-reversal-invariant 
momenta, the NM state obtained using the PBE-fSOC calcu¬ 
lation is found to be a QSH state. 

It is noticed that Pb atoms in the puckered TL lattice fa¬ 
vor the sp^ hybridization, leaving a single dangling bond 
on each Pb atom. Since the Dirac bands originating from 
such dangling-bond electrons have a small bandwidth of ^0.5 
eV Isee Fig. 2(a)], other electronic instabilities such as a 
charge/spin density wave (CDW/SDW) may be favored over 


the QSH state. For the CDW instability, we are unable to find 
its stabilization within the employed unit cell Isee Fig. 1(a)], 
but the ferromagnetic (FM) and AFM states can be stabilized. 
The calculated PBE and PBE-fSOC total energies of the FM 
and AFM states relative to the NM (or QSH) state are given 
in Table 1. We find that, for the PBE (PBE-fSOC) calculation, 
the FM state is less stable than the NM state by 6.8 (7.6) meV, 
while the AFM state is significantly more stable than the NM 
state by 37.4 (36.5) meV. 

TABLE I: Calculated total energies (in meV per unit eell) of the FM 
and AFM states relative to the NM state, obtained using the PBE and 
PBE+SOC calculations. The calculated spin magnetic moments (in 
Pb) of Pb and C atoms are also given. Here, three C atoms bonding 
to Pb are taken. 




PBE 

PBE+SOC 

FM 

AEpM-NM 

6.8 

7.6 


Mpb 

0.35 

0.33 


Me 

0.20 

0.21 

AFM 

AEafm-nm 

-yiA 

-36.5 


Mpb 

±0.32 

±0.31 


Me 

±0.16 

±0.22 


Figure 2(b) shows the PBE-fSOC band structure for the FM 
state. We find that the SOC opens a gap of Eg = 1.^ meV at 
the two points near the F point, while the PBE band structure 
shows a gap closing of the Dirac cones at Ef (see Fig. 5S(b) 
of the Supplemental Material IH]). From the PBE-fSOC cal¬ 
culation, we calculate the Berry curvatures O whose integral 
over the Brillouin zone (see Eig. 6S of the Supplemental Ma¬ 
terial IH]) gives the Chern number C = ^/bz^ 
with equal contributions from the valley regions. The nonzero 
Chern number is the characteristic of a QAH state, similar to 
the recently studied triphenyl-manganese lattice il. 

It is remarkable that AEM order is stabilized in a 
organometallic compound composed of only main group el¬ 
ements. Contrasting with the QSH and QAH states in the TL 
lattice, the AEM state exhibits a large band-gap opening at val¬ 
leys K and K' Isee Eig. 2(c)]: i.e., the PBE calculation gives 
an identical value of Eg = 469 meV, while PBE-i-SOC gives 
different values of 481 and 457 meV, respectively. This large 
gap opening of the AEM state gives a large thermodynam¬ 
ics stability relative to the NM state (see Table I). To under¬ 
stand the underlying mechanism for the AEM spin ordering, 
we plot in Eig. 2(d) the PBE spin-polarized local density of 
states (DOS) projected onto the Pb atoms at the A and B sub¬ 
lattices, together with their spin characters. It is seen that the 
occupied (unoccupied) spin-up and spin-down dangling-bond 
states are localized at the A(B) and B(A) sites, respectively. 
Since electronic states with the same spin direction can hy¬ 
bridize with each other, the hybridization takes place between 
the occupied and unoccupied spin-up (spin-down) states local¬ 
ized at the A(B) and B(A) sites, respectively, yielding not only 
a gain of the exchange kinetic energy but also an insulating 
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FIG. 2: (Color online) Band structures of the (a) NM, (b) FM, and (c) AFM states, obtained using the PBE+SOC calculation. The energy 
zero represents Ep. The insets in (a) and (b) magnify the energy dispersion near Ep. For the AFM state, the PBE spin-polarized local DOS 
projected onto the two Pb atoms at the A and B sites are given in (d). The charge characters of the spin-up (spin-down) states for the highest 
occupied and the lowest unoccupied bands are taken at the K point with an isosurface of 0.002 (—0.002) electrons/A^ 


gap Q. This kind of exchange interaction mediated by two 
electronic states which are energetically separated well below 
and above the Fermi level is characterized as a superexchange 
mechanism Issl [s^. It is likely that such a superexchange 
interaction can be facilitated due to a considerable hybridiza¬ 
tion of the Pb 6pz orbitals with the TT orbitals of neighboring 
phenyl rings. As shown in Fig. 7S of the Supplemental Mate¬ 
rial 1 I 33 I 1 and Table I, this hybridization is well represented by 
a large spin delocalization with the alternating spin densities 
of ^±0.32 and ~±0.16 pp for Pb atoms and C atoms in the 
phenyl ring, respectively. 


In Fig. 2(c), the SOC-induced renormalization of the val¬ 
ley gaps in the AFM state can be described by the parameter 
5: that is, = 2(m + t5), where m = ^Eg = 234.5 meV 
is the mass term arising from the AFM spin ordering and x 
= 1 (—1) represents the valley index of K {K'). The present 
value of 5 = 6 meV in the AFM state indicates a moderate 
spin-valley coupling (SVC), where the product v-x of spin and 
valley indices becomes a new degree of freedom of electrons. 
Here, the spin-up and spin-down bands remain degenerate be¬ 
cause the Hamiltonian is invariant under simultaneous time 
reversal T and spatial inversion P, although neither T nor P 
alone commutes with the Hamiltonian. Figure 3(a) (Fig. 8S 
of the Supplemental material ll^ ) shows the Berry curva¬ 
tures of the lowest (highest) unoccupied (occupied) spin-up 
and spin-down bands, computed using PBE-fSOC. It is seen 
that the spin-up and spin-down bands have the opposite values 
of ^2 at valley K and these values are also opposite to the cor¬ 
responding values at K', manifesting the SVC. Consequently, 
the Chern number becomes zero, and the insulating AFM state 
is characterized to be topologically trivial. 


Due to the SVC in the AFM state, the valleys’ degeneracy 
can be lifted without any external field |2ll[37]. For instance, 
the valleys can be doped asymmetrically with electrons (or 
holes), as schematically illustrated in Fig. 3(a). In the case 
of electron doping, the conduction band is occupied to shift 


(a) 


(b) 





(1/2,1) (-1/2,-1) 


FIG. 3: (Color online) (a) Schematic band structure of the AFM state 
with electron doping at valley K' and Berry curvatures of the lowest 
unoccupied spin-up and spin-down bands. The spin-up and spin- 
down currents under an applied in-plane electric field Ein are shown 
in (b), depending on the amount of doping level that shifts AEp < 
m + Ih and AEp > m - 1 - 25. Solid and dashed lines stand for the 
currents from K' and K, respectively. The spin and valley indices 
(s,x) are indicated in parentheses. 


Ep io 2 i higher energy. If this shift A^Ep is smaller than m 
+ 25, an applied in-plane electric field produce a transversal 
valley current with the electrons of (v,x) = (±^,—1), giving 
rise to a net transversal spin current without charge current, as 
schematically shown in Fig. 3(b). This is characterized as the 
QVH effect. It is notable that such spin currents are induced 
by the spin-dependent Berry curvatures Elj-t [see Fig. 3(a)] 
that determine the anomalous velocities of Bloch electrons, 
Na ^ Ein X For ^Ep > m + 25, the pair of currents 
with (v,x) = (±^,1) starts to contribute to a net transversal spin 
current in the reversal direction [see Fig. 3(b)]. When we 
switch from electron-doping to hole-doping, spin holes can 
generate valley and spin currents in the opposite directions. 
This generation of spin currents in the AFM state with time- 
reversal symmetry breaking represent the anomalous spin Hall 
effect. 

The mean field of the present AFM state can be char- 
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acterized as a spin-dependent staggered lattice potential in 
the honeycomb lattice. It is evident that an additional spin- 
independent staggered lattice potential will induce a closure 
of the gap in one of the valleys, eventually giving rise to a non¬ 
trivial electronic topology. To explore this topological phase 
transition, we apply an out-of-plane electric field Eout to pro¬ 
duce a staggered electric potential in the puckered honeycomb 
TL lattice. Here, Eout is simulated by superimposing an addi¬ 
tional sawtooth potential along the z direction with disconti¬ 
nuity at the mid-plane of the vacuum region of the supercell. 
This staggered electric potential is expected to lift the spin 
degeneracy of the AFM state because the occupied spin-up 
and spin-down DB states are localized at the A and B sites 
with different heights [see Fig. 2(d)], respectively. Conse¬ 
quently, as the magnitude of Eout increases, the gap of the 
spin-split bands would decrease, close, and reopen at the K 
and K' points, leading to a topological phase transition. In¬ 
deed, our PBE-fSOC calculations show the evolution of the 
band structure with increasing Eout, giving rise to a gap clos¬ 
ing of the spin-down band at valley K' for the critical field of 
Eout ~ 0-62 V/A. Figure 4(a) shows the band structure and 
Berry curvature obtained at Eout = 0-6 V/A. It is seen that the 
band gap at K {K') is much reduced as 119 (97) meV, com¬ 
pared to the value of 481 (457) meV in the absence of Eout 
[see Fig. 2(c)]. The calculated spin magnetic moments for 
the Pb atoms at the A and B sites become Mpi3(A) =0.28 and 
Mpb(B) = —0.25 }jiB, respectively, indicating a ferrimagnetic 
spin ordering. We note that the Chern number for this ferri¬ 
magnetic state is still zero. By contrast, for Eout = 0.7 V/A, 
the Chern number is changed into 1 as a consequence of the 
band inversion at valley K' where the sign of Berry curvature 
changes [see Fig. 4(b)]. Here, the spin magnetic moments for 
the Pb atoms are much reduced as Mpb(A) =0.18 and Mpb(B) = 
—0.15 }jiB. These results obviously indicate a phase transition 
to a topologically nontrivial QAH state with ferrimagnetic or¬ 
der. It is noticeable that, when Eout is above ^0.71 V/A, the 
gap of the spin-down band at valley K closes and reopens, 
which makes the system back to a trivial insulator. Since 45 
[i.e., the gap difference between K and K' valleys in Fig. 2(c)] 
amounts to 24 meV, the topological nontrivial QAH state in¬ 
duced by a certain out-of-plane external electric field may be 
realized at temperatures below ^300 K. 

In summary, we have performed comprehensive DFT cal¬ 
culations for the hexagonal TL lattice to investigate its true 
ground state by considering the QSH and QAH states as well 
as a band insulator with AFM order. We find that the topo¬ 
logically trivial AFM band insulator is more stable than the 
QSH and QAH states. This AFM state gives an emergent 
electronic degree of freedom characterized by the product 
of spin and valley indices, which leads to Berry curvature- 
induced spin and valley currents under electron or hole dop¬ 
ing. Furthermore, we find that the application of an out-of- 
plane electric field can induce a quantum phase transition from 
the AFM state to a topologically nontrivial QAH state with 
ferrimagnetic order. Our findings are fairly generic for 2D 
organometallic lattices, which will provide novel platforms 



FIG. 4: (Color online) Calculated band structures of (a) the topologi¬ 
cally trivial ferrimagnetic state at Eout = 0.6 V/A and (b) the topolog¬ 
ically nontrivial QAH state at Eout = 0.7 V/A. The Berry curvature 
of each structure is also given. 


for the applications of spintronics and valleytronics using 2D 
AFM hexagonal lattices. In particular, the metal centers and 
ligand choice in 2D organometallic lattices are two variables 
in their syntheses, which will surely enrich the exploration of 
topological and correlated electronic phases. 

This work was supported in part by National Research 
Foundation of Korea (NRF) grant funded by the Korean 
Government (2015R1A2A2A01003248), the National Nat¬ 
ural Science Foundation of China (Grants No. 11174009 
and No. 61434002), and the National Key Basic Research 
Program of China (Grants No. 2013CB921900 and No. 
2014CB921103). The calculations were performed by KISTI 
supercomputing center through the strategic support program 
(KSC-2014-C3-049) for the supercomputing application re¬ 
search. 

* Corresponding author: chojh@hanyang.ac.kr 
Corresponding author: jfengl 1 @pku.edu.cn 


[1] K. V. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45, 
494 (1980). 

[2] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005). 

[3] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005). 

[4] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Eirsov, Sci¬ 
ence 306, 666 (2004). 

[5] S. Murakami, Phys. Rev. Lett. 97, 236805 (2006). 

[6] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science. 314, 
1757 (2006). 

[7] M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, L. 
W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science. 318, 766 
(2007). 

[8] A. Roth, C. Briine, H. Buhmann, L. W. Molenkamp, J. Ma- 
ciejko, X.-L. Qi, and S.-C. Zhang, Science. 325, 294 (2009). 

[9] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Pang, and S.-C. 
Zhang, Nat. Phys. 5, 438 (2009). 


























5 


[10] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, 
D. Grauer, Y S. Hor, R. J. Cava, and M. Z. Hasan, Nat. Phys. 
5, 398 (2009). 

[11] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010) 
and references therein. 

[12] S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, 
S. von Molnar, M. L. Roukes, A. Y. Chtchelkanova, and D. M. 
Treger, Science 294 , 1488 (2001). 

[13] I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 323 
(2004). 

[14] A. Rycerz, J. Tworzydlo, and C. W. Reenakker, Nat. Phys. 3, 
172 (2007). 

[15] Y. P. Shkolnikov, E. P. De Poortere, E. Tutuc, and M. Shayegan, 
Phys. Rev. Lett. 89 , 226805 (2002). 

[16] D. Xiao, W. Yao, and Q. Niu, Phys. Rev. Lett. 99, 236809 
(2007). 

[17] X. Xu, Y. Wang, D. Xiao, and T. F. Heinz, Nat. Phys. 10, 343 
(2013). 

[18] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 1. 
Katsnelson, 1. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, 
Nature 438 , 197 (2005). 

[19] D. Pesin and A. H. MacDonald, Nat. Mater. 11, 409 (2012). 

[20] W. Yao, S. A. Yang, and Q. Niu, Phys. Rev. Lett. 102 , 096801 
(2009). 

[21] X. Li, T. Cao, Q. Niu, J. Shi, and J. Feng, Proc. Natl. Acad. Sci. 
110 , 3738 (2012). 

[22] D. Gunlycke, and C. T. White, Phys. Rev. Lett. 106, 136806 

( 2011 ). 

[23] Y. Jiang, T. Low, K. Chang, M. 1. Katsnelson, and F. Guinea, 
Phys. Rev. Lett. 110 , 046601 (2013). 

[24] Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, and M. 
S. Strano, Nat. Nanotech. 7, 699 (2012). 

[25] H. Zeng, J. DAi, W. Yao, D. Xiao, and X. Cui, Nat. Nanotech. 
7 , 490 (2012). 

[26] M. Ezawa, Phys. Rev. Lett. 110 , 026603 (2013). 

[27] M. Kim, C. H. Kim, H.-S. Kim, and J. Ihm, Proc. Natl. Acad. 
Sci. 109 , 671 (2011). 

[28] X. Qian, J. Liu, L. Fu, and J. Li, Science 346, 1344 (2014). 

[29] Z. F. Wang, Z. Liu, and F. Liu, Nat. Comm. 4 , 1471 (2013). 

[30] C.-C. Liu, W. Feng, and Y. Yao, Phys. Rev. Lett. 107 , 076802 
( 2011 ). 

[31] Z. F. Wang, Z. Liu, and F. Liu, Phys. Rev. Lett. 110, 196801 
(2013). 

[32] The present DFT calculations were performed using the VASP 
[Kresse and Hafner, Phys. Rev. B 48 , 13115 (1993); Kresse 
and Furthmiiller, Comput. Mater. Sci. 6, 15 (1996)] with the 
projector-augmented wave method and a plane wave basis set. 
For the treatment of exchange-correlation energy, we employed 
the generalized-gradient approximation functional of Perdew- 
Burke-Ernzerhof (PBE) [Perdew, Burke, and Ernzerhof, Phys. 
Rev. Lett. 77 , 3865 (1996)]. The TL lattice was modeled by a 
periodic slab geometry with ^25 A of vacuum in between the 
slabs. Here, we used the optimized lattice constant gq = 11.866 
A obtained using the PBE calculation with SOC. A plane wave 
basis was employed with a kinetic energy cutoff of 400 eV, and 
the k- space integration was done with the 9x9 Monkhorst-Pack 
meshes in the Brillouin zone. We employed a dipole correction 
that cancels the artificial electric field across the slab [Neuge- 
bauer and Scheffier, Phys. Rev. B 272 , 16067 (1992)]. 

[33] See Supplemental Material at 

http://link.aps.org/supplemental/xxxx for the PBE band 
structures of the NM, EM, and AEM states, the Berry curvature 
of the EM state, and the spin density of the AEM state. 

[34] K. Sato, L. Bergqvist, J. Kudrnovsky, P. H. Dederichs, O. Eriks¬ 


son, 1. Turek, B. Sanyal,G. Bouzerar,H. Katayama-Yoshida, V. 
A. Dinh, T. Eukushima, H. Kizaki, and R. Zeller, Rev. Mod. 
Phys. 82, 1633 (2010). 

[35] J. B. Goodenough, Phys. Rev. 100, 564 (1955). 

[36] J. Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 

[37] J. Qi, X. Li, Q. Niu and J. Peng, Phys. Rev. B 92, 121403(R) 
(2015). 


Supplemental Material for “Competing Magnetic Orderings and Tunable Topological 
States in Two-Dimensional Hexagonal Organometallic Lattices” 


Hyun-Jung Chaokai Ji Feng^’^^, Jun-Hyung Cho^*, and Zhenyu Zhang^ 

^ Department of Physics and Research Institute for Natural Sciences, Hanyang University, 17 Haengdang-Dong, 

Seongdong-Ku, Seoul 133-791, Korea 

^International Center for Quantum Design of Functional Materials (ICQD), Hefei National Laboratory for Physical Sciences at 
the Microscale (HFNL), and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of 

Science and Technology of China, Hefei, Anhui 230026, China 
^International Center for Quantum Materials, Peking University, Beijing 100871, China 
"^Collaborative Innovation Center of Quantum Matter, Beijing 100871, China 


1. Band structures of the NM, FM, and AFM states 

(a) (b) (c) 





MK' r KMMK' r KMMK F K'M 

FIG. 5S: Calculated band structures of the (a) NM, (b) FM, and (c) AFM states, obtained using the PBE ealculation. The inset in (a) magnifies 
the energy dispersion of the Dirac cone at the K point. In (b), the solid and dashed lines represent the spin-up and spin-down bands, respectively. 
The energy zero represents the Fermi level. 
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2. Berry curvature of the FM state 



FIG. 6S: (a) Calculated Berry curvature for the FM state, obtained using the PBE+SOC calculation. The Berry curvatures of the highest 
occupied state and the second highest occupied state are also given in (b) and (c), respectively. We find that the Chern numbers calculated from 
(a), (b), and (c) are 1, 0, and 1, respectively. 
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3. Spin density of the AFM state 



FIG. 7S: Tops and side views of the PBE spin density of the AFM state. The spin-up (spin-down) density is displayed in red (blue) eolor with 
an isosurface of 0.002 (—0.002) ^/A^. 
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4. Berry curvature of the AFM state 
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FIG. 8S: Calculated Berry eurvatures of the highest oecupied (a) spin-up and (b) spin-down bands in the AFM state, obtained using the 
PBE-i-SOC calculation. 
















